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ABSTRACT 

Dwarf galaxies tend to have redder [3.6 /im] — [4.5 /im] Spitzer broadband colors than spirals. 
To investigate this effect, for a large sample of dwarf galaxies we combine Spitzer fluxes with data 
at other wavelengths and compare to population synthesis models. Lower metallicity systems are 
found to have redder [3.6] — [4.5] colors on average, but with considerable scatter. The observed 
range in [3.6] — [4.5] color is too large to be accounted for solely by variations in stellar colors 
due to age or metallicity differences; interstellar effects must contribute as well. For the reddest 
systems, the 4.5 /im luminosity may not be a good tracer of stellar mass. We identify three 
factors that redden this color in dwarfs. First, in some systems, strong Bra emission contributes 
significantly to the 4.5 /im emission. Second, in some cases high optical depths lead to strong 
reddening of the starlight in the Spitzer bands. Third, in some galaxies, the nebular continuum 
dominates the 4.5 /im flux, and in extreme cases, the 3.6 /im flux as well. The harder UV radiation 
fields in lower metallicity systems produce both more gaseous continuum in the infrared and more 
Bra per star formation rate. The combination of these three factors can account for the 4.5 /im 
excess in our sample galaxies, thus it is not necessary to invoke a major contribution from hot 
dust to the 4.5 /im band. However, given the uncertainties, we are not able to completely rule 
out hot dust emission at 4.5 /im. More spectroscopic observations in the 3 — 5 /im range are 
needed to disentangle these effects. 
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The mid-infrared spectra of low metallicity 
dwarf galaxies are well-known to differ from those 
of spirals, with weaker polycyclic aromatic hy- 
drocarbon (PAH) e mission features and higher 
[Ne III] /[Ne ll] ratios (TThuan. Sauvage. fc Madder " 

~3lJS 



1999: 



Madden 2000; Galli ano et alj200 3: Houc k et al. 
Hunt. Bianchi. fe Maiolinol2005HMadden et all 
20061 ). This may be due to 



2004 

2006i IWu et al 
harder interstellar radiation fields (ISRFs) in 
dwarfs because of lower d ust extinction, lead- 
ing to PAH destruction ([Galliano et alj l2003l 



(|Galliano. Dwek. fe Chaiiiail l2008h . These dif- 
ferences are reflected in the mid-infrared broad- 
band colors of dwarfs, which show depressed 8 
/im e mission relative to 24 /im, compared to spi- 
rals dBoselli et all Il998t lEngelbracht et all [200l 
iRosenberg et al.ll2006l I2008T ). This is interpreted 
as a weakening of the 7.7 /im PAH feature in 
dwarfs. Dwarf and spiral spectral energy dis- 
tributions also differ at longer wavelengths; spi- 
ral spectra peak in the far-infrared, while dwarf 
spectra peak at shorter w a velengths, » 25 



60 urn jGalliano et al.1 120031: iHouck et all 12004 : 



Hunt. Bianchi. fc Maiolind 120051 ). This implies 



20051 : iMadden et al] l2006h . or to delayed ere 



ation of carbon dust by AGB stars in dwarfs 
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warmer large dust grains in dwarf galaxies on av- 
erage. Low metallicity galaxies also tend to have 
lower total dust-to-gas mass ratios than higher 
metallicity systems (e.g., Engelbracht et al. 2008; 



Hirashita et~aI1l2008fh 

The Spitzer [3.6] — [4.5] broadband colors 
of irrcgular/Sm galaxies also differ from those 
of sp i rals, being redder on average (IPahre et al 



2004; lEngelbracht et alj|2005t I Smith et al.ll2007[) . 



This color difference, however, has proven more 
difficult to interpret. Since the 3 — 5 /im ra- 
diation from galaxies may be a combination 
of starlight and interstellar emission, it can be 
difficult to accurately disentangle the contribu- 
tions from these components to th e observed 
broadband fluxes. IPahre et al.l ( 2004 ) suggested 
that this reddening is due to younger stars in 
frr/Sm galaxies compared to spirals, with hot- 
ter stars having redder [3.6] — [4.5] colors. In 
general, hotter stars tend to have redder L — 
M co lors (~[3.6] — [4-5]) than cooler stars (|Cox 
2000). Alternatively, these color differences have 
been attributed to hot dust in dwarf galax- 



ies c o ntributing at 4.5 /xm (jEngelbracht et al 



20051: iHunter et all 120061 ). The Bra line may 



also contribute significa n tly to the 4.5 urn ban d 
(|Churchwell et al.l l200i lElmegreen et al.1 120061 ). 
In addition, the nebular continuum in the in- 
frared can also contribute significantly at wave- 
lengths > 3 ixm in low metallicity systems (e.g., 
Kriiger. Fritze-v. Alvensleben. fc Lo ose 1995). An- 



other factor may be the 3.3 /zm PAH feature, 
which contributes to the 3.6 fim Spitzer band. 
A lower heavy element abundance, and therefore 
less PAHs, can redden the [3.6] — [4.5] color. Also, 
stars with low metal abundances may have redder 
colors because of less line blanketing in the 4.5 /im 
band (M. Cohen 2006, private communication). 
Evolved stars can also redden this color, for ex- 
ample, Mir as are f — 2 magnitudes redder in L — 
M th an non- variable optically-bright stars (jSmith 



2003). In highly obscured systems, reddening due 



to dust extinction can also be important e ven in 
the Spitzer bands (e.g., Roussel et al.ll2006l ). 

To determine the origin of the red [3.6] — [4.5] 
colors in dwarf galaxies, for a large sample of 
dwarfs we have collected Spitzer mid-infrared 
fluxes, metallicities, and ground-based optical, 
near-infrared, and radio measurements from the 
literature. Combining these data with population 
synthesis modeling, we have investigated trends in 
the stellar and interstellar spectra with metallicity 
and age. In Section 2 of this paper, we describe 
the sample and the data. In Section 3, we compare 



the [3.6] — [4.5] colors with other properties of the 
galaxies. We discuss expected [3.6] — [4.5] colors 
from starlight in Section 4, and in Section 5, we 
investigate Bra contributions to the 4.5 /im band. 
In Section 6, we discuss reddening due to dust 
extinction in the Spitzer bands. In Section 7, we 
discuss reddening of the [3.6] — [4.5] color due to 
the nebular continuum. In Section 8, non-stellar 
contributions to the 3.6 /im Spitzer band are in- 
vestigated. In Section 9, the question of hot dust 
contributions at 4.5 /im is discussed. Morphology 
is discussed in Section 10, while luminosity and 
stellar mass are investigated in Section 11. Some 
additional issues are discussed in Section 12, and 
conclusions are presented in Section 13. 

2. The Sample and the Data 

To investigate how the [3.6] — [4.5] color de- 
pends upon metallicity and other properties of 
dwarf galaxies, we collected Spitzer 3.6, 4.5, 5.8, 
8.0, and 24 /xm fluxes from the literature for a 
sample of 71 irregular, blue compact dwarf (BCD), 
Sm, Sd, and other dwarf galaxies, for which oxy- 
gen abundances are available. This sample is not 
complete in any sense, however, it spans a large 
range of metallicity, absolute magnitude, Spitzer 
luminosities, and gas content. Table 1 lists the 
references from which we extracted the data used 
in this study. 

For all of the galaxies in our sample with very 
red or very blue [3.6] — [4.5] published colors, 
we downloaded the Spitzer archival images and 
checked the published values. For all except the 
bluest galaxy in the sample, UGC 4483, the col- 
ors we obtained agreed with the published colors 
within ^0.15 magnitudes. For UGC 4483, how- 
ever, our [3.6] — [4.5] color is —0.16, compared 
to -0.39 from lEngelbracht et al.l (|2008t ). For this 
galaxy, we used our fluxes rather than the pub- 
lished numbers. 

We also collected total-galaxy 21 cm HI and 
Ha fluxes for the sample galaxies. We also 
used Ha/H/3 ratios from the literature, preferably 
global (total-galaxy) values when possible; other- 
wise, we used smaller aperture measurements. In 
addition, we searched for both total-galaxy and 
smaller aperture Bra and Br7 measurements. We 
also obtained H/3 equivalent widths and measure- 
ments of the 9.7 /im silicate absorption feature 
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optical depth for our sample galaxies, preferring 
global values over small aperture measurements 
when available. In Table 1, for each type of data, 
we give the number of galaxies with such data 
available. 

We also obtained total-galaxy optical UBVRI 
and near-infrared JHK broadband fluxes for our 
sample from the literature. These magnitudes 
were corrected for Galactic extinction, when nec- 
es sary, and converted to Jy usin g the zero points 
in Bessel. Castelli. fe Pled (1998). In a number of 



cases, large discrepancies were found between dif- 
ferent published magnitudes for the same galaxy. 
In these cases, we did not use either measurement 
in our analysis. 

As a comparison sample of spiral galaxies, 
we started with the sample of nearby 'normal' 
(not strong l y inte racting) spirals discussed in 



(not strong l y inte racting) spirals discussed m 
ISmith et al.l ( 2007 ). These galaxies have total- 



galaxy Spitzer infrared and ground-based op- 
tical and near-infrare d broa dband fluxes avail- 
( 20071) . Most h ave total 



able from Dale et al 



Ha fluxes available from iKennicutt et al.l (|2003l ). 



and Pact /Ha ratios for th e inner 50" available 



from ICalzetti et al.l (|2007l ). As spirals tend to 
have radial abundance gradients, for compar- 
ison to the dwarfs we selected spirals which 
have disk-averaged oxygen abundances avail- 
able ( Moustakas fe Kennicuttll2006t Calzetti et al 



2007). After removing Sd galaxies, there are 13 



galaxies in the spiral sample. 

3. [3.6] [4.5] Colors as a Function of 
Metallicity, Gas Content, and Star For- 
mation Rate 

In the left panel of Figure 1, we compare the 
[3.6] — [4.5] color against oxygen abundance for 
our sample dwarf galaxies. The right panel shows 
the same plot for the spirals. The dwarfs are on 
average redder in [3.6] — [4.5], in addition to hav- 
ing lower metallicities. For the dwarfs, a relation is 
seen, in that colors redder than [3.6] — [4.5] > 0.3 
are only seen at low abundances, log(0/H) + 12 < 
8.2 (< 0.5 solar). In contrast, the spirals all have 
[3.6] - [4.5] - 0.0, and disk-averaged log(0/H) + 
12 > 8.2. 

For the dwarfs, there is considerable scatter in 
Figure 1. To test whether this scatter is due solely 
to differing amounts of interstellar matter, with 



more gas-rich galaxies having a larger amount of 
dust and therefore more 4.5 /^m excess, in Figure 1 
we color-coded the data points according to their 
m(HI) — [3.6] color, where m(HI) is the magni- 



tudc in the 21 cm HI line (de Vaucouleurs et al 



Il99lt IPaturel et aT1l2003l ). The corresponding HI 
mass for the dwarfs range from 7 x 10 6 M Q to 7 
x 10 9 Mq. In this plot, as log(0/H) increases, 
m(HI) — [3.6] tends to increase (i.e., the HI mass 
to 3.6 fim luminosity ratio M#j/L 3 . 6 decreases). 
This is consistent with res ults of previous studies 
(e.g.. lFisher fc Tullvlll975h . which show that lower 
metallicity systems tend to be richer in HI relative 
to their stellar mass. However, for a given metal- 
licity, no strong trend of m(HI) — [3.6] with [3.6] 
— [4.5] is found. Thus the 4.5 fim excess does not 
appear to be due to an increasing mass fraction of 
interstellar matter. This result is uncertain, how- 
ever, because molecular gas is not included in this 
plot. Many of these galaxies are undetected in CO, 
the standard tracer of molecular gas in galaxies, 
and the amount of molecular gas present is uncer- 
tain. 

In Figure 2, we compare the [3.6] — [4.5] col- 
ors of these galaxies with the [3.6] — [24] colors. 
The left panel shows the dwarfs, the right the spi- 
rals. In this Figure, we have color-coded the data 
points as a function of oxygen abundance. The 
24 /jm luminosity is considered a good tracer of 
star formation rate, as it arises from small dust 



grain s heated by UV photons (e.g.. ICalzetti et al 



2005). If the 3.6 /mi flux is mainly due to older 
stars, then [3.6] — [24] is roughly a measure of the 
mass-normalized star formation rate, and there- 
fore the intensity of the UV field. Figure 2 shows 
that indeed more intense UV fields are correlated 
with more 4.5 /Ltm excess. The distributions of the 
spirals and the dwarfs overlap in Figure 2, but the 
dwarf sample extends to much redder [3.6] — [24] 
colors as well as redder [3.6] — [4.5] colors than the 
spirals. The dwarf galaxies with the redder colors 
tend to have lower metallicities, however, some low 
metallicity galaxies have blue Spitzer colors. Thus 
low metallicity appears to be a necessary but not 
sufficient condition for red [3.6] — [4.5] and [3.6] 
- [24] colors. 
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4. Stellar [3.6] [4.5] Colors from Popu- 
lation Synthesis 

To test whether stars alone can account for 
the variations in [3.6] — [4.5] colors, we used 
the Starburst99 stellar population synthesis code 
" 1999h to produce model stel- 
these bands. We used version 
which includes the Padova 



(jLeitherer et al 



lar colors in 
5.1 of this code 
asymptotic giant branch ste llar evolution models 
(| Vazquez fc Leithererl 120051 ). We note that this 
code does not model red supergiants in low metal- 
licity systems well ([Vazquez fe Leitherer 2005h . 
which introduces some uncertainty to our results. 
We used a grid of models with ages that range 
from 1 Myrs to 20 Gyrs, and metallicities of 1/50, 
1/5, 1/2.5, 1, and 2.5 times solar. 

Dwarf galaxies have a large range of mass- 
normalized star formation rates; some are under- 
going violent bursts of star formation, wh ile oth- 
ers have low current rates (jHunterl Il997l ). The 
integrated UBV colors of dwarf irregular galax- 
ies are general ly consistent with constant star for- 
matio n rates (jHunter fc Gallagher! 11985c Ivan Zee 
20011 ). Blue compact dwarf galaxies are undergo- 
ing recent starbursts, but an older population is 
usual l y also present dShulte-Ladbeck et al. 1998, 



20011 iDrozdovskv et all 120021 : iKong et alJ 12003 : 
Thuan fe Izotovll2005f ). As limiting cases to the 
true star formation histories of the galaxies, we 
therefore ran both continuous and instantaneous 
burst star formation models. 

We integrated the Starburst99 model spectra 
over the Spitzer bandpasses and calculated the 
[3.6] — [4.5] colors. In Figure 3, we plot model 
[3.6] — [4.5] colors vs. age for various metallicities. 
In Figure 3, we only include stellar contributions 
to the Spitzer fluxes; no dust emission or interstel- 
lar gas emission lines are included. The models 
shown in Figure 3 also do not have reddening due 
to dust extinction included (see Section 6). Color 
corrections are also not included, as these are ex- 
pected to be very small (see IRAC Data Manual). 
In the top panel, we present instantaneous burst 
models, while the bottom panel shows continuous 
star formation models. 

Figure 3 shows that, according to Starburst99, 
there is a relatively small spread in the [3.6] — 
[4.5] colors of stars of different ages and metal- 
licities, <0.2 magnitudes. As expected, there is 



a trend such that younger (<6 x 10 6 Myrs) and 
lower metallicity systems (except >2.5 Z Q ) have 
redder colors than intermediate age systems (<10 8 
Myrs). However, the spread in colors is quite 
small, and cannot account for the full range of 
observed colors for dwarfs. Except for very young 
(<5 Myrs) models with greater than solar metal- 
licity, the model stellar colors are between —0.16 
and 0.04. Of our dwarf galaxies, 70% are redder 
than this limit. Thus the observed range in [3.6] — 
[4.5] colors of dwarfs is not due solely to variations 
in star colors due to age and metallicity variations. 

5. Bra Contributions to the 4.5 fim Spitzer 
Band 

In this section, we estimate the importance 
of the Bra line at 4.05 /jm, which may con- 
tribu te significantly to the Spitzer 4.5 /im band 
(e.g.. lChurchwell et al.ll2004l ). For example, in the 
Galactic H II region M17, Bra is estimated to con- 
tribute 20% of the total 4.5 /jm broadband flux 
( Povich et al.l 120071 ) . At low metallicities, stars 



are hotter and therefore more ionizing photons 
are available for the same mass of stars. Thus as 
metallicity decreases, for the same star formation 
rate the Ha and Bra fluxes are expected to in- 
creas e (e.g., ICharlot fc Longhetti 2001 : Lee et al 
2002f ). 



To determine how much Bra can affect the 
[3.6] — [4.5] color, we added Bra to the Star- 
burst99 models shown in Figure 3. We used the 
Ha fluxes provided by Starburst99, along with a 
no minal Bra/Ha ratio of .0273 (10,000 K, Case 
B: lHummer fc Storevlll987l ). We then added this 
model Bra flux to the Starburst99 stellar spec- 
trum and integrated over the Spitzer bandpasses 
to get model colors. In Figure 4, we plot the [3.6] — 
[4.5] model color including Bra against age, with 
different metallicity models plotted with different 
symbols. These models do not include dust ex- 
tinction or emission, or other interstellar lines be- 
sides Bra. A comparison of Figure 4 to Figure 3 
shows that, in some very young systems, Bra can 
contribute significantly particularly for low metal- 
licity systems. For example, for ages < 10 Myrs 
and a metallicity of 1/50 solar, Bra may redden 
the [3.6] — [4.5] color by >1 magnitude. A higher 
assumed temperature will decrease this reddening 
somewhat. 
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A more direct way to determine the Bra con- 
tribution to the observed 4.5 /im flux is to ei- 
ther measure it directly, or to extrapolate from 
B1-7 or Ha measurements. Only one galaxy in 
our sample, II Zw 40, had a total-galaxy Bra 
flux available (jVerma et al.ll2003l ) . From this mea- 
surement, we find that only ~11% of the Spitzer 
4.5 /zm flux of II Zw 40 is due to Bra. Another 
galaxy in our sample has a total-galax y published 
Bry flux available (SBSG 0335-052; I Dale et al " 



2001). Using the smal l- aperture spectros copic 



pect 

Bra/Br7 flux ratio from | Vanzi et al.l ([2000) and 
Hunt. Vanzi. fe Thuanl ( 2001 ) to extrapolate to 



the total Bra flux, we estimate that the Bra con- 
tribution to the 4.5 /im Spitzer flux for this galaxy 
is ~31%. This may be an over-estimate, since 
some of the Br7 flux may aris e in less extincted 
regions (e.g.. lVanzi et alJl2000h . 

For NGC 5253, the B r a me asurement of 
Kawara. Nishida. fe Phillips! (| 19891 ) in a 10" x 



20" aperture gives a lower limit to the Bra con- 
tribution to the Spitzer 4.5 //m flux of 2%. This 
aperture is considerably smaller than the optical 
angular size for NGC 5253 of ~2' x 5' given in 
the NASA Extragalactic Database (NED) and the 
observable extent of the 3.6 /im emission in the 
archival Spitzer image of ~1' x 2'. However, the 
infrared emission is strongly centrally peaked in 
the archival Spitzer images, with 25% of the total 
3.6 /im flux and 47% of the 4.5 /jm flux within a 
10" x 10" region. Thus the Bra contribution to 
the 4.5 /im flux is likely not much more than 4 — 
8% for NGC 5253. 

For galaxies with published total-galaxy Ha 
measurements, we estimate the Bra contribution 
indirectly from these fluxes, correcting for inter- 
nal extinction using the published Ha/H/3 ratio. 
From the extinction-corrected Ha flux, we calcu- 
lated the Bra flux, and determined the fraction of 
the observed 4.5 /im broadband flux due to Bra. 
This calculation shows that, if the Ha/H/3 esti- 
mates of extinction are appropriate, for our sam- 
ple galaxies the fraction of the observed Spitzer 
flux due to Bra ranges from 0.3% to 23%, with 
most (80%) having less than 10% of their total 4.5 
/zm flux due to Bra. In some cases, Ha/H/3 may 
underestimate the extinction (see Section 6), thus 
these percentages are lower limits. 

We then re-calculated the [3.6] — [4.5] colors of 
the galaxies, after removing the flux due to Bra. 



In the left panel in Figure 5, for the dwarf galaxies 
with published Ha fluxes, we plot the metallicity 
against both the observed and the Bra-removed 
[3.6] — [4.5] colors. For the galaxies with measured 
Bra or Br7, we plot two Bra-corrected values, one 
calculated using the Ha (filled red squares) and 
one using the Bra or Br7 value (open green trian- 
gles). In Figure 5, we mark the range in [3.6] — 
[4.5] expected from starlight alone, —0.16 < [3.6] 
— [4.5] < 0.04 (see Figure 3). The correction for 
Bra shifts the galaxies to the left, to bluer [3.6] — 
[4.5] colors, with the shifts ranging from — 0.3 
magnitudes. Note that for II Zw 40 and SBSG 
0335—052, the Bra correction inferred from Ha 
and Ha/H/3 measurements is less than that deter- 
mined from a direct Bra measurement (in the case 
of II Zw 40) or from Br 7 (for SBSG 0335-052). 
This suggests that the Ha/H/3 ratio may under- 
estimate the extinction in some cases (see Section 
6). 

Most of the shifts seen in Figure 5 are relatively 
small, <0.1 magnitudes. However, for a few of our 
galaxies (for example, He 2-10), the shifts are suffi- 
cient to move the galaxies into the range expected 
from starlight. In these cases, Bra may be mainly 
responsible for the red [3.6] — [4.5] colors. 

However, for some of our sample galaxies (e.g., 
II Zw 40, SBSG 0335-052, HS 0822+3542, NGC 
4194, and SHOC 391), the Bra-corrected colors 
are still redder than those expected by starlight 
alone. For II Zw 40, this conclusion is particularly 
robu st, because a direct Bra measurement is avail- 
able (IVerma et al. l l2003l). For SBSG 0335-052, the 
Br7 observation of lDale et alJ (|200lh also provides 
a reasonably strong constraint. 

In the right panel of Figure 5, a similar plot 
is shown for the spiral sample. The open black 
circles are the published [3.6] — [4.5] values, while 
the red crosses show the values after correction 
for Bra, calculated using the total-galaxy Ha 
fluxes and the Pao/Ha ratio for the inner 50" 



(jCalzetti et all 120071) . This plot shows that for 



normal spirals the corrections for Bra in the 
Spitzer 4.5 /im filter is generally very small, con- 
sistent with the fact that most of the published 
[3.6] — [4.5] values are in the expected range for 
starlight. 
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6. Reddening of Starlight in the Spitzer 
Bands due to Dust Extinction 

Another possible source of reddening of the [3.6] 
— [4.5] color is dust extinction of the starlight. 
The model [3.6] — [4.5] colors of stars shown in 
Figure 3 do not include reddening due to dust ex- 
tinction. In modeling the spectral energy distribu- 
tions in galaxies, reddening is sometimes neglecte d 
at wavelengths > 3.6 /zm (e.g. JPraine et al.ll2007l h 
However, if the infrared sources are highly ob- 
scured, reddening by extinction may be important, 
even in the Spitzer bands. 

How much reddening is present depends upon 
the properties of the dust grains, and therefore 
the extinction law, as well as the geometry of the 
system. The ratio of visual extinction to B — V 
reddening Ay/E(B — V) = Ry has been found to 
vary from ~ 3.1 in diffuse gas to ~5.5 in dens e 
clouds (e.g., ICardelli. Clayton, fc Mathisl 1 19891 ). 
The extinction law i n the Spitzer bands als o varies 



with optical depth ( Chapman et al. 2008?) . Com- 



bining the extinction laws from these papers gives 
a range for Ay/E(3.6 — 4.5) between ~ 41 for 
diffuse gas to ^68 for more extincted regions. 

The published Ha/H/3 ratios for our sample 
dwarfs imply Ay between and 2.5, or redden- 
ing in [3.6] — [4.5] of less than 0.06 magnitudes, 
not sufficient to account for the 4.5 /jm excess in 
our reddest systems. However, there is evidence 
that the Ha/H/3 ratio strongly underestimates the 
extinction in at least some of our systems. 

The Ha/H/3 ratio of SBSG 0335-052, the galaxy 
in our sample with the reddest [ 3.6] — [4.5] 
color , implies only Ay ~ 0.3 — 0.5 (llzotov et al 



1997t llzotov. Thuan. fc Stansihskal l2007f ). How- 
ever, mid-infrared spectroscopy of SBSG 0335-052 
shows a deep 9.7 iim silicate abso r ption feature 
( Thuan. Sauvage. fc Maddenl Il999l : iHouck et al" 
20041 ). suggesting a much higher optical depth of 
Ay ~ 15 — 20. The Br7/Bra ratio also suggests 
a much larger extinction t han that implied by the 
Ha/H/3 ratio, Ay ~ 12 (jHunt. Vanzi. fc Thuan 
200 lh . Dust modeling of the 6 - 100 /jm 
spectral energy distribution also suggests that 
most of the mid-infrared emission originates 
from a highly obs cured (Ay ~ 12—30) huni- 
nous star cluster (Hunt. Vanzi. fc Thuanl 12001 



Hunt. Bianchi. fc Maiolinql2005l : IPlante fc Sauvaga 



rived extinction with wavelength implies a more 
complex geometry than a simple uniform fore- 
ground screen of dust, and suggests that the 
light seen in the optical is less obscured than the 
bulk of the starburst 



Hunt. Vanzi. fc Thuan 



20011 : iReines. Johnson, fc Huntl 120081) . The 3-4 



/im so urce in SBSG 0335-052 is ve ry compact, 
<1.2" (jHunt. Vanzi. fc Thuanl l200ll ). consistent 
with this picture. If the extinction obtained from 
the silicate absorption, Ay ~ 15—20, applies to the 
majority of the stars contributing to the 3.6 and 
4.5 /im emission, then E(3.6 - 4.5) ~ 0.22-0.48. 
Combined with Bra, this can account for some, 
but not all, of the observed reddening in [3.6] — 
[4.5]. 

In the left panel of Figure 6, we plot the oxygen 
abundance of the sample dwarf galaxies against 
the [3.6] — [4.5] color, after this color has been cor- 
rected for both Bra and reddening of the starlight 
due to dust extinction. For these galaxies, we plot 
extinction corrections made using the 9.7 /jm sili- 
cate optical depth (blue crosses), the Bra/Bi-7 ra- 
tio (magenta asterisks), and/or the Ha/H/3 ratio 
(black filled diamonds), depending upon the avail- 
able data. As in Figure 5, we have marked the 
expected colors due to starlight alone. Note that 
the [3.6] - [4.5] color of SBSG 0335-052 is too red 
to be accounted for by stars, even after correcting 
for both Bra and extinction. Thus there must be 
another contributor to the 4.5 /jm excess for this 
galaxy. 

A similar situation may exist f or II Zw 40 . 
The Ha/H/3 ratio implies A y ~ 1.1 (|Frenchlll98oT : 
Kinman fc Davidson! Il98l! k or E(3.6 - 4.5) ~ 



0.03. However, longer wavelength data gives 
highe r extinctions. The Bra/Br/3 ratio gives Ay 
~ 10 ( Verma et al.ll2003l) . while comparison of Br7 
to radio continuum implies Ay ^8 — 10 and mod- 
eling of the infrared s pectral energy distribution 
gives Ay ~ 20 - 30 (jHunt. Bianchi. fc Maiolinol 
2005! ). As with SBSG 0335-052, the mid-infrared 



source in II Zw 40 is very compact, with most 
emission confined to 0.5" (Beck et al. 2002). In 
contrast to SBSG 0335-052, however, no silicate 
absorption feature is seen in mid-infrared spec- 



tra of II Zw 4 ([Martin-Hernandez et all l2006t: 



Wu et al. 20061). thus the extinction is uncertain. 



200a iTakeuchi et all 12003). The increase in de- 



If Ay ~ 10 holds for the bulk of the population 
contributing to the Spitzer broadband flux, then 
E(3.6 - 4.5) - 0.15 to 0.24. As shown in Figure 
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6, this could account for some, but not all, of the 
observed reddening. 

NGC 5253 is another example of a galaxy 
with a larger implied optical depth from the 
Bra/Br7 ratio (Ay ~ 11 in a 10 " x 20" beam; 
Kawara. Nishida. fe Phillips! Il989h and the sili- 
cate absorption (A y ^ 8 — 25 in a 5" 4 aperture; 



galaxies different datasets are available. 

In the right panel of Figure 6, the blue open 
diamonds show the [3.6] — [4.5] colors of the nor- 
mal spirals after correction for both Bra and for 
starlight reddening. The latter te rm was calc u- 



Aitken et al. .1982) than from the Ha/H/3 ratio 



(Ay ~ 0.3 for the whole galaxy; iMoustakas fc Kennicuttl 
20061 ). As with II Zw 40 and SBSG 0335-052, ex- 



tinction and Bra together cannot account for all 
of the reddening in [3.6] — [4.5] (see Figure 6). 
As noted in Section 5, however, our correction for 
Bra may be slightly too small. 

For He 2-10 (ESO 495-21), if the silicate absorp- 
tion ( Ay - 15 in a 5'/9 aperture; IPhillips. Aitkens. fe Rochfre a bove range in Ay for spirals (see also|Calzetti 



lated using the Paa/Ha ratios of ICalzetti et al 
((2007J). Note that these corrections are small 
for the spirals, with Ay ranging from 0.2 — 2.5. 
Thus the correction for starlight reddening, like 
that for Bra, does not shift the spirals much 
on this plot. We note that other techniques for 
obtaining extinctions such as background counts 
of galaxies through galaxian disks also give rela- 
tively small average obscurations for nearby spirals 
( White. Keel, fc Conselicell2000l ). consist ent with 



1984} or the Bra /Br7 -implied extinction (Ay ~ 17 



in 7" 1 x 3" 5 aperture: iKawara! Nishida. fc Phillips! 
1989) are used to correct the Spitzer colors, in- 



2001). 



stead of the exti nction implied by th e Ha/H/3 
ratio (Ay ~ 1.7; IVacca fc Contil Il992h . the cor- 
rected [3.6] — [4.5] color becomes too blue (see 
Figure 6). In this galaxy, the lower extinction 
estimate may be more appropriate. 

Two of our other dwarfs with red [3.6] — [4.5] 
colors, NGC 4194 and Haro 11 (ESO 350-38), also 
show apparent sil icate absorption in t heir mid- 
infrar ed spectra ( Brandl et al.1 l2006t IWu et al 
20061 ) . For Haro 1 1 , no silicate optical depth anal- 



ysis has been done yet. For NGC 4194, the sili- 
cate optical depth i s estimated to be r g .7 = 0.37 
Brandl et all 120061). Assuming Ay/A 9 . 7 - 18.5 
Roche fc Aitkenlll984h gives E(3.6 - 4.5) between 



0.10 and 0.17. Combined with the correction for 
Bra, this can account for most, but not all, of the 
excess 4.5 ^im emission (see Figure 6). 

In contrast to th ese galaxies, I Zw 18 shows no 
silicate absorption ( Wu et aD 12007 ) and the ex- 
tinction implied by t he Ha/H/3 ratio is very small 
( Cannon et al. 20021 ). In this case, the correction 
for extinction is small compared to that for Bra 
(see Figures 5 and 6). 

One caveat in this analysis is that the redden- 
ing of the stars in a galaxy may differ significantly 
from that of the ionized gas, depending upon the 
geometry of the system. Thus caution should be 
taken in applying extinction corrections obtained 
from emission line ratios or silicate absorption fea- 
tures to stellar colors. Furthermore, for different 



7. Contributions from the Nebular Con- 
tinuum 

Another source of near- to mid-infrared light 
in galaxies is nebular continuum emission, which 
we have neglected so far in our analysis. In the 
case of low metallicity star burst systems, this 

component can be important dJov fc Lester 19881 

Kriiger. Fritze-v. Alvensleben. fc Loosell995 : Vanzi et al 



200oilHunt. Vanzi. fc Thuanll200l!) . To dete 



the contribution from the gaseous continuum to 
the Spitzer bands, we ran another set of Star- 
burst99 models, including just starlight and the 
gaseous nebular continuum, excluding Bra and 
dust reddening. To calculate the gaseous contin- 
uum, Starburst99 uses parameters from Ferlandl 
(1980). The longest wavelengths used in these 
tables is 4.5 fim, therefore the predicted values 
of the nebula r continuum are unre liable at longer 
wavelengths ( Leitherer et al. I ll999l) . 

In Figure 7, we plot the model [3.6] — [4.5] 
colors as a function of age for various metallici- 
ties. Comparison with Figure 3 shows that the 
addition of nebular continuum emission signifi- 
cantly reddens the [3.6] — [4.5] colors of young sys- 
tems, particular for low metallicities. For example, 
<10 Myr 0.02 Z instantaneous burst models with 
the nebular continuum included have [3.6] — [4.5] 
colors about 0.5 magnitudes redder than models 
without this component. For solar metallicities, 
ages <4 Myrs produce such red colors. Thus the 
nebular continuum may be a very important com- 
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ponent in very young and/or low metallicity sys- 
tems. 

In extreme cases, the nebular continuum may 
dominate not only the 4.5 fim. band, but also 
the 3.6 /urn band. This is illustrated in Figures 
8 and 9, where we plot the ratio of the nebu- 
lar continuum flux in the 3.6 /urn and 4.5 /xm 
bands to the flux from starlight plus the neb- 
ular continuum, for the models shown in Fig- 
ure 7. Figures 8 and 9 show that, for solar 
metallicity and ages < 3 Myrs, and for 0.02Z Q 
and < 10 Myrs, contributions to both bands 
from the nebular continuum are larger than from 
starlig ht. iKriiger. Fritze-v. Alvensleben. fc Loose! 
(|l995l ) found that, for a 5 Myr O.O2Z burst, -70% 
of the flux in the L band is due to the nebular 
continuum, consistent with these plots. Since the 
nebular contributions are typically larger at 4.5 
/im than 3.6 /im, the net effect is a reddening of 
the observed colors of galaxies. The issue of non- 
stcllar emissions at 3.6 fim is discussed further in 
Section 8. 

To determine the contribution from the nebular 
continuum for individual galaxies in our sample, 
we require an estimate of the average age of the 
stellar population. We do this two ways. First, we 
use Starburst99 models and published H/3 equiva- 
lent widths to estimate the age. For each galaxy, 
we use models with a metallicity closest to that of 
the published metallicity, and we fit both instan- 
taneous and continuous star formation models. 
Second, when sufficient broadband optical data is 
available (>3 filters), we fit the optical spectral en- 
ergy distribution to Starburst99 models to get the 
age and extinction of the contributi ng stars, and 



the uncertainties on these v alues fe.g. JSmith et al 



20081 : lHancock et all 120081 ). We include the nebu- 
lar continuum in the optical in doing these fits, and 
we conservatively assumed uncertainties in the op- 
tical colors of 0.1 magnitudes for all galaxies. In 
convolving the model spectra with the response 
functions of the broadb and filters, we adde d the 
Ha emission line (e.g., ISmith et al.l 120081 ). but 
did not include other emission lines. We used 
the lCardelli. Clavton. fc Mathid (|l989h extinction 
law. 

We did not use the near-infrared photometry 
in the fits, but used them after the fact to decide 
whether instantaneous or continuous star forma- 
tion is more appropriate for the system. In most 



cases, both instantaneous and continuous models 
were consistent with the near-infrared data; in a 
few galaxies, the continuous models produced too 
much near-infrared or 3.6 /im flux. 

We then used the best-fit ages from the instan- 
taneous burst models to determine the amount 
of reddening in the [3.6] — [4.5] color due to the 
nebular continuum. We then corrected the ob- 
served colors by this amount, adding to the cor- 
rections for Bra and dust extinction discussed 
above. As discussed in Section 4, some dwarf 
galaxies can be fit by continuous star formation 
models, while others are better explained by mul- 
tiple bursts. For normal spirals, the broadband 
optical spectral energy distributions can gener- 
ally be fit by co ntinuous or exponentially decaying 
star formation dSearle. Sargent, fc Bagnudol 19731 : 
lLarson fc Tinslevill978l ). We choose to use instan- 
taneous burst models to calculate the nebular con- 
tribution to the 4.5 /im band since they produce 
less nebular emission for the same age, as shown 
in Figures 8 and 9. This thus provides a limiting 
case for the nebular continuum contribution. For 
some galaxies, the fit to the optical data was quite 
poor, likely because their star formation history 
does not match either an instantaneous or con- 
tinuous burst. These galaxies were excluded from 
Figure 10. 

In Figure 10, the nebular-corrected [3.6] — [4.5] 
colors are plotted against oxygen abundance, with 
the dwarfs in the left panel and the spirals in 
the right. In these plots, the filled upside-down 
magenta triangles are the nebular continuum- 
corrected colors calculated using the best-fit ages 
from the H/3 equivalent widths, while the cyan 
open diamonds are the nebular-corrected values 
using the best-fit ages from the broadband opti- 
cal fitting. To illustrate the uncertainties in this 
correction, the open red circle shows the color af- 
ter correction for nebular emission using the best 
fit age from the broadband fitting minus the la 
uncertainty in the age. In some cases, the correc- 
tions to the Spitzer [3.6] — [4.5] color due to the 
nebular continuum is quite large, for example, it 
is ~0.4 magnitudes for SBSG 0335-052. 

In most cases, within the uncertainties the com- 
bined contributions of Bra, reddening due to dust 
extinction, and the nebular continuum can ac- 
count for all of the observed reddening above 
starlight. In Figure 10, for all of the galaxies ex- 
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cept two, the two estimates of the corrected color 
(the magenta upside-down filled triangle and the 
cyan open diamond) or the lower limit on the cor- 
rected color (open red circle) lie within the green 
hatched region, or span this region. For the two 
exceptions, the datasets are incomplete, thus we 
cannot make firm conclusions. For Haro 11, as 
noted in Section 6, silicate ab sorption may b e 
present in the Spitzer spectrum ( Wu et al. 2006h . 
however, no optical depth has been published, and 
no published Bra or Br7 values are available. For 
SHOC 391, no Bra, Br7, or silicate absorption 
measurements are available. Thus the extinction 
of the starlight and/or the correction for Bra in 
the Spitzer 4.5 /im band may be under-estimated 
for these galaxies. 

8. Non-Stellar Contributions to the 3.6 
/im Band 

It is often assumed that the 3.6 /im band is 
dominated by starlight, and a standard [3.6] — 
[4.5] spectrum is assumed for the stellar compo- 
nent, which is then subtracted from the observed 
4.5 1 MB, flux to get the r esidual 4.5 /xm emission 



(e.g.. iHelou et alJ l2004t lEngelbracht et all 12005k 



Dale et alJl2005l) . However, as noted above, nebu- 



lar continuum emission may be important at both 
3.6 /im and 4.5 /im (see Section 7). Also, if hot 
dust is significant at 4.5 /zm, it ma y also be presen t 
at 3.6 /im, and even at 2 /im (e.g.. lLu et al.l 2003). 

To test whether significant non-stellar emis- 
sion is present at 3.6 /jm in individual galaxies 
in our sample, we used our Starburst99 fits to 
the broadband optical data. For all but one of 
the galaxies with available optical data, we can- 
not rule out that all of the 3.6 /zm emission is due 
to starlight, since our best-fit stellar spectra are 
consistent with the observed Spitzer 3.6 /im flux 
within the uncertainties. The only galaxy in our 
sample that shows a clear 3.6 /jm excess above the 
inferred stellar continuum is the reddest galaxy in 
[3.6 ] - [4.5], SBSG 0335-052 (s ee Figure 11; also 
see iGalliano. Dwek. fc Chaniall l2008i For SBSG 
0335-052, although the inferred stellar continuum 
is consistent with the near-infrared photometry for 
both the instantaneous and continuous bursts, it 
is about a factor of 5 below the Spitzer 3.6 /zm 
measurement. However, when the nebular contin- 
uum is included in the Spitzer bands, most of the 



3.6 /im flux can be accounted for (see Figure 11). 
This shows that the 3.6 /im excess may be mainly 
due to the gaseous continuum. 

We note that the best-fit extinction for the 
Starburst99 model of SBSG 0335-052, A v = 1.0 
± 0.5, is significantly less than that inferred by the 

9.7 (jm silicate absorption fe a ture ( Ay ~ 15 — 20; 



Thu an. Sauvage. fc Madden 19991 : Houck et al 



2004). As mentioned in Section 6, for this galaxy, 
stars that dominate the observed light in 



the 



the optical may be much less obscured than the 
stars that dominate in the infrared. Thus our 
Starburst99-derived extinction based on the opti- 
cal data may not be an appropriate value to use 
in the infrared. A second more highly obscured 
stellar population, unseen in the optical, may also 
contribute to the observed Spitzer 3.6 /im flux of 
this system. However, this second population is 
not evident in the near-infrared data. 

Unfortunately, for some galaxies in our sample, 
no suitable broadband optical data were available 
for fitting, or we were not able to get good fits to 
the available data with our Starburst99 models. 
Thus we were not able to use this method to test 
for 3.6 /zm excesses above the stellar continuum 
in these galaxies. These galaxies include some of 
the reddest galaxies in our sample in [3.6] — [4.5]: 
SHOC 391, Haro 11, Pox 4, and HS 0822+3542. 
Further observations and more detailed modeling 
of these systems are needed to search for 3.6 /zm 
excesses above the stellar continuum. 

None of the spirals in our sample for which a 
good fit to the optical data was found showed a 3.6 
/xm excess above the inferred stellar continuum. 
As they also have [3.6] — [4.5] colors consistent 
with starlight (Figure 1), this also implies little 
4.5 /im non-stellar excess in the spirals. As noted 
earlier, the 3.3 /im PAH feature lies within the 3.6 
/im Spitzer band, however, the spirals in our sam- 
ple do not show unusually blue [3.6] — [4.5] colors. 
Thus there is no strong evidence for a strong PAH 
contribution to the 3.6 /im flux of the spirals in our 
sample. As noted earlier, dwarf galaxies tend to 
have weaker PAH features than spirals, thus this 
conclusion also may hold for the dwarf galaxies. 

An alternative way to search for 3.6 /im excesses 
is to compare the [3.6] — [4.5] colors with other 
near-to- mid-infrared colors (Figure 12). The spi- 
rals (right panel) show less scatter in these colors 
than the dwarfs (left panel) , and the higher metal- 
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licity dwarfs (blue open squares) show less scatter 
than the lower metallicity dwarfs. For the dwarfs, 
weak correlations are seen between [3.6] — [4.5] 
and K — [3.6], and between [3.6] — [4.5] and H — 
[3.6], in that the extremely red galaxies in [3.6] — 
[4.5] also tend to be red in these other colors. This 
suggests that non-stellar contributions at 3.6 /xm 
may be present in these systems. The correlation 
seen between [3.6] — [4.5] and K — [4.5] is not 
unexpected, if K roughly traces the older stellar 
population. No correlation is seen between [3.6] — 
[4.5] with H - K. 

9. Contributions at 4.5 /im from Hot In- 
terstellar Dust 

It is often assumed that excess emission at 4.5 



Lim i s due to hot dust (e.g., lEngelbracht et al 



20051 iHunter et aD l2006h . However, Figure 10 



shows that, for all the galaxies in our sample with 
sufficient data available, a combination of Bra, 
the nebular continuum, and dust reddening may 
be able to account for most of the 4.5 /zm excess. 
Thus a major contribution from dust emission is 
not required. However, with the available data, 
given the uncertainties we are not able to com- 
pletely rule out hot dust emission in the 4.5 /zm 
band. 

One way to distinguish between these differ- 
ent processes is with spectroscopy. However, at 
the present time, few high S/N spectra of galaxies 
in the 3 — 5 /zm range are available, since the 
Spitzer spectrometers do not reach wavelengths 
< 4.9 /im. Co-added Infrared Space Observatory 
(ISO) 2.5 — 4.9 /zm spectra of disk galaxies show 
a continu um excess wit h a color temperature of 
~1000K (iLu et al.ll2003[h At longer wavelengths, 
in the ~24 /im range, the continuum emission from 
star forming g alaxies is attributed to 'very sma ll 
grains' (VSGs: lDesert. Boulanger. fc Pugedll990h . 
In Galactic HII regions, the 6 — 16 /urn VSG con- 
tinuum is strong, in contrast to the PAH features 
which weaken in HII region s ( Cesarskv et aljl996l 
IPeeters. Spoon, fc Tiele'nl \2004 ). Within dwarf 
galaxies, the ratio of the PAH feature strength to 
the VSG continuum varies spatially, being anti- 
correlated with the [Ne III 15.8 /zm]/[Ne II 12.8 
/im] line ratio (jMadden et al. 2006; Beirao et al 



to vary much within a dwarf galaxy, this indicates 
that PAHs are destroyed in hard UV fields, while 
mid-infrared VSG excitation is enhanced. 

In a few low metallicity galaxies with spatially- 
resolved ISO spectra, the short wavelength end 
(4.9 — 5.6 /im) shows strong conti nuum emission 
in reg ions with strong Ha emission (jMadden et al 



2006). This supports the idea that the VSG con- 
tinuum may extend well into the 3 — 5 /im range. 
For SBSG 0335-052, a Sp itzer spectrum show s a 5 
/xm continuum ~ 2 mJy ([Houck et al.ll2004T ). Ex- 
trapolation of this to shorter wavelengths could 
possibly account for much of the observed 4.5 
/zm broadband emission. Three of our galaxies 
with very red [3.6] - [4.5] colors, SBSG 0335- 
052, Haro 11, and SHOC 391, have been classi- 
fied as 'mid-infrared peakers' by lEngelbracht et al 



(2008), having infrared spectral energy distribu- 
tions that peak near 24 /im rather than in the 
far-infrared. This is consistent with these galaxies 
having very hot dust grains on average, and the 
dust continuum extending down to <4.5 /im. 

Thus with the available data we cannot rule out 
hot dust as a contributing factor. However, the 
above discussion shows that the nebular contin- 
uum can also be quite strong at these wavelengths, 
thus it is hard to separate these two components 
in the 3 — 5 /im wavelength range. 

10. Morphology 

Of the 71 galaxies in our dwarf sample, 18 
are classified in the NASA Extragalactic Database 
(NEE0) as blue compact dwarfs or 'compact', 19 
as irregular (Im/10/Irr/IBm/IB/IAB), 3 as Sm, 
and 3 as Sd, with the rest listed as other types, 
such as HII galaxies or starbursts. In some cases, 
the classifications in NED differ from those in the 
SIMBAD database!!, as dwarf classifications are 
often subjective, and these catalogs are heteroge- 
neous. 

In the left panel of Figure 13, we plot the 
observed [3.6] — [4.5] colors of the dwarfs vs. 
log(0/H) + 12 as in Figure 1, however, in Fig- 



2006). Since [Ne III]/[Ne II] is a measure of radia- 



tion field hardness, and metallicity is not expected 



lr The NASA/IPAC Extragalactic Database (NED) is oper- 
ated by the Jet Propulsion Laboratory, California Institute 
of Technology, under contract with the National Aeronau- 
tics and Space Administration. 

2 http:simbad.u-strasbg.frsimbad, operated at CDS, Stras- 
bourg, France 
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ure 13 the points are color-coded according to the 
morphological type in NED. The blue filled cir- 
cles are classified as blue compact galaxies or as 
compact, while the magenta open diamonds are 
listed as irregular. The green crosses are galax- 
ies with other classifications. In the right panel 
of Figure 13, we plot [3.6] - [24] vs. [3.6] - [4.5] 
for the dwarfs as in Figure 2, however, the points 
are color-coded according to morphological type 
rather than metallicity. 

Figure 13 shows that the galaxies with the red- 
dest [3.6] — [4.5] and [3.6] — [24] colors are more 
likely to be classified as blue compact dwarfs than 
as irregular. However, there is considerable scat- 
ter in this plot, with some compact galaxies hav- 
ing blue Spitzer colors. The observation that the 
reddest galaxies tend to be classified as compact 
agrees with the hypothesis that more compact, 
and therefore more UV- intense and/or more ob- 
scured, star formation is occurring in these sys- 
tems. 



1989), and interpreted as a mass-abundance rela- 
tion. For irregulars, this relation has been ex- 



tended to the 4.5 fjxa luminosity by iLee et af 



( 20061 ). who found less scatter at 4.5 /xm than at 
B, and concluded that the 4.5 /xm band is a better 
tracer of stellar mass than the B band. 

In the left and right panels of Figure 14, we 
compare the 4.5 /mi luminosity with log(0/H) + 
12 for our sample dwarfs and normal spirals, re- 
spectively. For consistency with previous work, 
we use the 'monochromatic' luminosity jaL„, with 
frequency 6.67 x 10 13 Hz. The luminosities were 
calculated using distances from NED calculated 
with H Q = 73 km s _1 Mpc -1 , and Virgo/Great 
Attractor/Shapley cluster infall, except in the 
case of NGC 6822, where this method fails. For 
NGC 6822, we use the Cepheid distance from 



Paturel et al. ( 20021 ). As expected, the 4.5 /xm lu 



The reddest systems also s how evidence for in- 
terac tions or mergers: I Zw 18 (Ivan Zee. Skillman. fc Salze^ NED . as in Figure 13 
19981 ). SB SG 0335-052 dPustilnik et all l200ll). 



minosities of the dwarf sample are typically much 
lower in than for the spirals, but there is some 
overlap in range. The points in Figure 14a are 
color-coded according to their morphological type 



Haro 11 JOstlin et alj l200ll). HS 082 2+3542 
dPustilnik et al.l l2003t ICorbin et all 120061) . NGC 



In Figure 14, we include the best-fit line from 



Lee et al.l (|2006l ) for their irregular galaxies. Fig- 



4194 ( Armus, Hcckma n. fc Milevlll990D. and II Zw 
40 (|Baldwin. Spinrad. fc Terlevichlll982h . This is 
consistent with the idea that an interaction or 
merger has driven gas into the center of the sys- 
tem, causing very centrally-concentrated star for- 
mation. This may lead to very intense UV fields 
and/or high optical depths, causing hotter dust 
grains, more 24 /xm emission, more reddening due 
to extinction, more gaseous continuum emission, 
and/or more Bra. 

In our earlier Sp itzer study of mor e massive in- 
teracting galaxies (|Smith et al.l 120071 ) we did not 
find a statistically significant difference between 
their [3.6] - [4.5] colors and the colors of normal 
spirals. However, that sample was selected to be 
early-stage encounters, with the disks widely sep- 
arated, while some of the dwarfs in the current 
sample may be more advanced interactions. 



ure 14 shows that most of our sa mple galaxi e s, and 
all of our irregulars, lie near the ILee et al. (2006) 
line. However, a handful of galaxies show large 
excesses in 4.5 /xm luminosity for their metallic- 
ities. These include some of the systems noted 
earlier for having very red [3.6] — [4.5] colors: 
SBSG 0335-052, SHOC 391, Harol 11, and I Zw 
18. In contrast, some of the galaxies with red 
[3.6] — [4.5] colors, su ch as II Zw 40, lie close 
to the ILee et al.l (|2006l ) line. In addition, a few 
other galaxies that are less extreme in [3.6] — [4.5] 
also show up as discrepant in Figure 14: SHOC 
567, Tol 2138-405, and UM 420. Tol 2138-405 
and SHOC 567, like the redder galaxies, have 
eculiar morphologies indicative of interactions 



Telles. Melnick. fc Terlevichlll997l : iKniazev et al 



20041 ) . UM 420 lies behind a foreground ellipt i- 



cal galaxy ( Salzer. MacAlpine. fc Boroson 1989), 
thus the published Spitzer flux may be confused. 

Figure 14 shows that for extreme starbursts, 
compact systems, and/or interacting systems, the 
4.5 /xm luminosity may not be a good tracer of stel- 
lar mass. This may be in part due to non-stellar 
contributions in the 4.5 /im band. Alternatively, 
abundance has been observed (e.g.JSkillman, Kennicutt, %Hqdj&action or merger may have driven unpro- 



11. Spitzer Luminosities and Stellar Masses 

For irregular and dwarf elliptical galaxies, a 
correlation between blue luminosity and oxygen 
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cessed gas into the central region of these galaxies, 
lowering t heir metallicity relative to their stellar 
mass (e.g-. lRupke. Veilleux. fe Bakerlliooih . 

In Figure 15, we plotted the 3.6 /im luminosity 
vs. oxygen abundance for both the dwarfs and the 
spirals. Again, a correlation is seen, in that higher 
luminosity systems tend to have higher metallic- 
ities. However, as in Figure 14, there are some 
galaxies with too high luminosities for their metal- 
licities. These are the same galaxies that are dis- 
crepant in Figure 14, thus these galaxies may also 
have large non-stellar excesses at 3.6 /im, or may 
have lowered metallicities due to interactions or 
mergers. 

12. Discussion 

For the dwarf galaxies in our sample, we find 
that lower metallicity systems tend to have redder 
[3.6] — [4.5] colors than higher metallicity systems, 
in that systems with log(0/H) + 12 < 8.2 tend to 
be redder on average. Comparison to stellar pop- 
ulation synthesis models shows that the observed 
colors of many dwarfs are too red to be accounted 
for solely by unobscured starlight. 

For all of the galaxies in our sample with suf- 
ficient data available, a combination of Bra emis- 
sion, reddening by dust, and nebular emission can 
account for the observed 4.5 /im excess. Nebular 
continuum emission, which is often neglected in 
discussions of the Spitzer colors of galaxies, may 
contribute to both the 3.6 /zm and 4.5 /zm bands, 
and in extreme cases can dominate the observed 
fluxes in these bands. As seen in Figure 2, 4.5 
/zm excess is correlated with the mass-normalized 
24 (im luminosity and therefore with more intense 
UV fields. In lower metallicity systems, the 4.5 
/zm excess is larger for the same [3.6] — [24] color, 
consistent with both stronger Bra emission and 
more nebular continuum for lower abundances. 

Another factor that is sometimes ignored is red- 
dening by dust extinction in the Spitzer bands, 
which can be quite large in some situations. Some 
of our sample galaxies show deep silicate absorp- 
tion features and/or large Bra/Br7 ratios, in spite 
of the small extinctions implied by Ha/H/3 or 
optical colors. In these cases, reddening of the 
starlight due to dust extinction may be quite im- 
portant. The situation in these very obscured 
dwarfs may resemble that in the starburst ellip- 



tical galaxy NGC 1377, which has a deep sili- 
cate absor ption feature and a high implied opti- 
cal depth (|Roussel et alJl2006t ). NGC 1377 may 
contain a very young highly obscured compact 
starburst ( Roussel et al.|[2006h . Like our extreme 
dwarfs, NGC 137 7 has a very red [3 6] - [4.5] color 
(see Figure 16 in lSmith et aDl2007t ). 

Given the uncertainties in the data and models, 
we cannot rule out an additional contribution to 
the 4.5 /jm flux from dust heated to high temper- 
atures by very intense UV fields. This dust con- 
tinuum may be the short-wavelength extension of 
the VSG continuum seen at longer wavelengths in 
ISO and Spitzer spectra of Galactic HII regions 
and star forming regions in dwarf galaxies. At the 
present time, the nature of the grains responsible 
for this 4.5 /zm dust emis sion is unclear. Ac cord- 
ing to the dust models of Draine & Lj| (l20f)7h . the 
grains responsible for the continuum at 6 — 30 
/im have sizes ~ 15— 40A , and may be large PAHs 
(N c ~ 2000 - 3 x 10 4 ). braine fc Lil (120071) con- 
cluded that in most cases this continuum emission 
is due to single-photon heating, however, in the 
case of very intense UV radiation fields, a parti- 
cle may not cool completely before absorbing an- 
other photon. In this case, higher temperatures 
are reached, and the peak of the VSG spectrum 
moves to shorter wavelengths. In sources with 
strong 4.5 /im continuum, the particles may be 
larger and/or hotter, and possibly may be in ther- 
mal equilibrium with a hard intense UV field. 

It is unclear how the grains that produce the 
strong 6 — 30 /im continuum in HII regions are 
related to the carrie rs of the ~2 — 3.7 /im con- 



tinuum found by ISellgren. Werner, fc Dinerstein 
([19831 ) in reflection nebulae, which was attributed 
to small grains (~10A) stochasti cally heated t o 
-T000K by UV or optical photons (|Sellgrenlll984h . 
The ambient UV fi elds in these tw o envir onments 
are quite different. An fc Sellgrenl ( 2003 ) suggest 
that the near-infrared continuum-producing par- 
ticles in reflection nebulae may be PAHs that are 
larger, more ionized, and/or more dehydrogenated 
than those that produce the 3.3 /im PAH feature, 
or may require higher energies for excitation. 

Interestingly, a 4.5 /im excess has also been 
found for the diffus e Galactic interst ellar medium 
Draine fc il (|2007t ) found that the iFlagev et a! 



([20061 ) Spitzer broadband measurements of diffuse 
Milky Way dust are consistent with standard dust 
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models, except for the 4.5 p, m flux, which is too 
high by a factor of ~_ 1.6. Flagev et al. ( 20061 ) 



and iDraine fe Lil (|20071 ) conclude that there is a 
'near- infrared continuum' whi ch contributes at 4.5 
/im and 3.6 /im. Furthermore. iLi fc Draind (|200lh 
compared their dust models to high Galactic lati- 
tude emission as seen by DIRBE, and found a 4.9 
/im excess (~2.5x). They suggested an additional 
form of opacity at these wavelengths, for example, 
carbon chains (e.g., Allamandola et al. 1999). It is 
unclear whether the particles that produce this 4 
— 5 (im excess in the diffuse ISM are the same as 
those that are responsible for the excess in dwarf 
galaxies. 

( 20051) 



As noted by lEngelbracht et al 



and 



Draine et all (|2(307[ )7 there appears to be a sharp 



change in the 8/24 /im flux ratio of galaxies at 
log(0/H) + 12 = 8.2 (-1/3 solar), in that this ra- 
tio abruptly decreases at this metallicity, perhaps 
because of weakening of the 7.7 /mi PAH feature 
at low metallicities. We see an abrupt change in 
the [3.6] — [4.5] color near this same metallicity, 
suggesting that the 4.5 /mi excess and the PAH 
deficiency are related. One possibility is that at 
this metallicity threshold the UV field hardens suf- 
ficiently to more efficiently destroy PAHs, and at 
the same time produces more Bra, more infrared 
nebular continuum, and/or more effectively excite 
the carriers of 4.5 /im dust emission. A second 
possibility is that the destruction of the carriers of 
the PAH emission features produces the particles 
that are responsible for a 4.5 /im dust continuum. 

A third possibility is that there is an abrupt 
change in the chemistry of the interstellar mat- 
ter at this metallicity, leading to both a defi- 
ciency in PAHs and to a larger percentage of grains 
that are more effective at radiating in the 4.5 /im 
regime. For ex ample, the galaxian chemical evo - 
lution model of lGalliano. Dwek. fc Chanial ( 20081 ) 
shows a time delay between the production of sil- 
icate dust grains in supernovae and the formation 
of carbon-rich material in asympotic giant branch 
stars, which can account for the relative deficiency 
of PAHs at low metallicities. As noted earlier, 
the galaxies in our sample with extremely red 
[3.6] — [4.5] colors have infrared spectral energy 
distributions that peak near 2 4 /im, while other 



low m etallicity systems do not. lEngelbracht et al 



in these galaxies, or are inefficiently heated. Thus 
the composition of the interstellar dust in these 
galaxies may be different from other low metallic- 
ity systems. 

For all of our galaxies with sufficient optical 
data except one, population synthesis modeling 
cannot rule out that most of the 3.6 /mi Spitzer 
flux is due to the stars seen at optical wavelengths. 
The exception, SBSG 0335-052, either has a highly 
obscured stellar population unseen at optical and 
near-infrared wavelengths, or a strong nebular or 
hot dust continuum that contributes significantly 
at 3.6 /im as well as at 4.5 /im. A similar situa- 
tion applies for a star forming clump in the tidal 
tail of the interacting galaxy pair Arp 285. For 
this clump, population synthesis based on optical 
broadband data indicates an excess at both 3.6 /im 
and 4.5 /im flux above the inferred stellar contin- 



( 20081 ) suggest that the particles responsible for 
far-infrared/submillimeter emission are deficient 



uum (|Smith et al.ll2008D . As with SBSG 0335-052, 
this may be due to either a strong nebular contin- 
uum, hot dust being present at 3.6 /im, or high 
obscuration in the optical. 

13. Conclusions 

For a large sample of dwarf galaxies, we show 
that the Spitzer broadband [3.6] — [4.5] colors de- 
pend upon metallicity, in that lower metallicity 
systems (log(0/H) — 12 < 8.2) show redder col- 
ors on average than higher metallicity systems, but 
with much scatter. Comparison with stellar popu- 
lation synthesis models shows that starlight alone 
cannot account for the range of observed colors, 
even when accounting for variations in age and 
metallicity. Many galaxies have an excess at 4.5 
/im compared to that expected from stars. In some 
galaxies, Bra contributes significantly to the ob- 
served 4.5 /im flux. In other systems, the nebular 
continuum plays an important role in reddening 
the [3.6] — [4.5] color, and in extreme cases may 
dominate the observed flux at both 3.6 /im and 
4.5 /xm. The lower metallicities in dwarfs hardens 
the UV field, increasing the Bra luminosity and 
the gaseous continuum in the infrared. In addi- 
tion, high obscuration can also redden the [3.6] — 
[4.5] color in some cases. The combined effects of 
these three factors can account for the observed 
[3.6] — [4.5] colors of all the galaxies in our sam- 
ple that have sufficient data available. Given the 
uncertainties, we cannot rule out an additional 
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component from hot dust at 4.5 fim. Harder UV 
fields may increase the average dust temperature, 
increasing dust contributions at 4.5 //m. However, 
for our sample galaxies it does not appear to be 
the dominant factor. 

The galaxies with the largest 4.5 /im excess 
are classified as blue compact dwarfs, and show 
evidence of recent interactions and mergers with 
other galaxies. Such encounters may produce 
very concentrated star formation, causing very in- 
tense UV fields and/or high extinctions. For these 
galaxies, the 4.5 //m luminosity may not be a good 
tracer of stellar mass. 

A similar analysis was done for a sample of 
nearby normal spirals. This showed that contri- 
butions from Bra and the nebular continuum to 
the 4.5 /im fluxes of normal spirals are generally 
small, as is reddening of the [3.6] — [4.5] starlight 
by dust. This is consistent with their observed 
[3.6] — [4.5] colors of ~ 0.0. For normal spirals, 
these bands appear to be dominated by starlight. 
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Fig. 1. — Left: A plot of log(0/H) + 12 vs. [3.6] — [4.5] for the sample dwarf galaxies. The different colors 
represent different HI mass to stellar mass ratios, as indicated by the HI magnitude m(HI) minus the 3.6 ^m 
magnitude [3.6]. The HI magnitude is defined as in Paturel et al. (2003). The magenta open diamonds are 
the most Hi-rich systems, with mag(HI) - [3.6] < 1.3. Red filled circles have 1.3 < mag(HI) - [3.6] < 3.0, green 
crosses 3.0 < mag(HI) - [3.6] < 4.1, cyan filled squares 4.1 < mag(HI) - [3.6] < 5.5, and blue asterisks 5.5 < 
mag(HI) - [3.6]. Galaxies marked in black open squares do not have HI data available. Typical uncertainties 
on the Spitzer [3.6] — [4.5] colors of galaxies are ^0.05 (Smith et al. 2007). Some of the galaxies are labeled. 
Right: A similar plot for the comparison sample of 'normal' spiral galaxies. 
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Fig. 2. — Left: The [3.6] — [24] vs. [3.6] — [4.5] colors for the sample dwarf galaxies. The different colors 
represent different metallicity ranges. The red open triangles are log(0/H) + 12 < 7.8, the green filled circles 
7.8 < log(0/H) + 12 < 8.2, and the blue open squares log(0/H) + 12 > 8.2. Right: A similar plot for the 
'normal spirals' sample. 
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Fig. 3. — Left: The [3.6] — [4.5] colors for the Starburst99 models as a function of age. Only starlight is 
included in these models; no contributions from the interstellar medium are included. The top panel shows 
the results for the instantaneous burst models; the lower panel gives the continuous star formation models. 
The different colors and symbols represent different metallicities, with red (filled squares) being 1/50 Z©, 
green (open circles) 1/5 Z©, blue (crosses) 1/2.5 Z©, cyan (open diamonds) 1 Z©, and black (filled circles) 
2.5 Z©. Right panel: lower two plots are the same as the left panel, except with an expanded x-axis, for 
comparison with Figures 1 and 2. The top panel on the right shows a histogram of the [3.6] — [4.5] colors 
for the dwarf sample, for comparison with the models. 




[3.6] - [4.5] 



Fig. 4. — Model [3.6] — [4.5] colors for stellar population models as a function of age, as in Figure 3, except 
with Bra added to the 4.5 /im band. The symbols and colors are as in Figure 3. Note that the scale on the 
x-axis extends further than in Figure 3b. 
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Fig. 5. — Left: The oxygen abundance of the subset of our sample dwarf galaxies with published Ha 
fluxes, plotted against the observed [3.6] — [4.5] color (open black squares) and the [3.6] — [4.5] color after 
correction for Bra. The filled red squares are the Bra-corrected colors after correction using Ha; the open 
green triangles are after correction using measured Bra or Br7 fluxes. The Bra-corrected datapoints are 
connected to the observed datapoints by dotted lines. The expected range for unextincted starlight is shown 
by hatch marks. Right: A similar plot for the sample spirals, where the open cir cles are the unco r rected 
points, and the red crosses are after correc tion for Bra us i ng to tal Ha fluxes from iKennicutt et al. (2003) 



and Paa /Ha ratios for the inner 50" from ICalzetti et al 
published Paa/Ha values are plotted. 



(2007). Only galaxies with total Ha fluxes and 
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Fig. 6. — Left: The same as Figure 5a, except with points added to show an additional correction for 
reddening of the starlight due to extinction by dust. Blue crosses represent points corrected using extinction 
estimates from the 9.7 /im silicate absorption feature. Magenta asterisks show points corrected using Bra/Br7 
ratios, and black filled diamonds show points corrected for reddening due to dust extinction using Ha/H/3 
ratios. The correction for extinction has been added to the correction for Bra in this plot. The extinction- 
corrected points are connected to the Bra-corrected points by dotted lines. Right: The same as the left 
panel, except for the spiral sample. The b lue open diamonds sh ow the colors after correction for reddening 
of the starlight using the Paa/Ha ratio of Calzetti et al. ( 20071 ). 
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Fig. 7. — Model [3.6] — [4.5] colors as a function of age as in Figure 3, except with the nebular continuum 
added. No Bra or reddening by dust extinction is included. The symbols are as in Figure 3. The x-axis 
range has been set to match that in Figures 1,2, and 3b. 
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Fig. 8. — The results of Starburst99 models showing the ratio of the flux in the 3.6 /im band from the nebular 
continuum to that from starlight plus the nebular continuum. The different colors and symbols represent 
different metallicities, with red (filled squares) being 1/50 Z , green (open circles) 1/5 Z , blue (crosses) 
1/2.5 Z , cyan (open diamonds) 1 Z , and black (filled circles) 2.5 Z . 
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Fig. 9. — The results of Starburst99 models showing the ratio of the flux in the 4.5 /J,m band from the nebular 
continuum to that from starlight plus nebular continuum. The symbols are as in Figures 3, 4, 7, and 8. No 
Bra is included. 
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Fig. 10. — As in Figure 6, except with points added to show an additional correction for reddening due to 
the nebular continuum. The left panel shows the results for the dwarfs, the right the spirals. The magenta 
filled upside-down triangles show the nebular continuum-corrected points, calculated using ages determined 
from the H/3 equivalent width. The cyan diamonds show the nebular continuum-corrected points calculated 
using the best-fit age from the broadband optical data. The red open circles shows the nebular-corrected 
values, calculated using the best fit age from the broadband data, minus the ler uncertainty to the age. The 
correction for the nebular continuum has been added to the corrections for Bra and extinction in this plot. 
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Fig. 11. — The optical-mid-infrared spectral energy distributions of SBSG 0335-052. The black solid curve 
is the best-fit Starburst99 model, while the blue dot-dashed curve shows the lower limit to the age, and 
the red dashed curve the upper limit to the age. The upper and lower limits curves are plotted with their 
best-fit extinctions. The green dotted curve shows the best-fit curve including the nebular continuum in the 
Spitzer bands. The other curves only include starlight in the Spitzer bands; for these curves, the nebular 
continuum was only included in the optical and near-infrared. The lower panel shows the instantaneous 
burst models, while the continuous models are shown in the upper panel. These models were calculated 
using a metallicity of l/50th solar. The models are scaled to the V band optical flux. The Ha and Bra lines 
are plotted for comparison, but no other emission lines . The optic al, near-infrared, and Sp itzer data are 
from lPustifnik. Pramskii. fc Kniazekl ( 2004 ). Vanzi et al. ( 2000l ). and Engelbracht et al. (2008), respectively. 
The best-fit ages are 37 ± Jjl Myrs for the continuous burst model, and 6 ± | Myrs for an instantaneous 
burst. 
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Fig. 12. — Left: Various near-to-mid-infrared color-color plots for our dwarf sample galaxies. The red open 
triangles are log(0/H) + 12 < 7.8, the green filled circles 7.8 < log(0/H) + 12 < 8.2, and the blue open 
squares log(0/H) + 12 > 8.2. Right: Same plots for the spiral sample. 
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Fig. 13. — Left panel: Same as Figure la, with the observed [3.6] — [4.5] colors for the dwarfs as a function 
of oxygen abundance, but with the data points color-coded according to morphological types from NED. 
The blue filled circles are classified as BCD galaxies or as compact, while the magenta open diamonds are 
listed as Irr, Im, 10, IB, or IAB. The green crosses are galaxies classified as other types such as HII galaxies 
or starbursts. Right panel: Same as Figure 2a, with the [3.6] — [4.5] color plotted against [3.6] — [24], except 
with the points color-coded according to NED morphological type instead of metallicity. 
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Fig. 14. — Left panel: the 4.5 ixm luminosity vXj v plotted against oxygen abundance log(0/H) + 12, for 
our sample dwarfs. The I Lee fc Skillmanl (|2004[ ) best-fit line for irregular galaxies is shown. Right panel: a 
similar figure for our spirals. The data points are color-coded according to morphological types from NED 
as in Figure 13. The blue filled circles are classified as BCD galaxies or as compact, while the magenta open 
diamonds are listed as Irr, Im, 10, IB, or IAB. The green crosses are galaxies classified as other types such 
as HII galaxies or starbursts. 
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Fig. 15. — Left panel: the 3.6 fxm luminosity v\j v plotted against oxygen abundance log(0/H) + 12, for 
our sample dwarfs. Right panel: a similar figure for our spirals. The data points are color-coded according 
to morphological types from NED as in Figure 13. The blue filled circles are classified as BCD galaxies or 
as compact, while the magenta open diamonds are listed as Irr, Im, 10, IB, or IAB. The green crosses are 
galaxies classified as other types such as HII galaxies or starbursts. 
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Table 1 

Sources of Data used in this Study 



Type of Data 



Number 
Available 



References 



Spitzer 3.6, 4.5, 5.8, 8.0, and 24 Jmx fluxes 7T 
Oxygen Abundances 71 



21 cm HI fluxes 



Ha fluxes 



40 



40 



Ha/H/3 ratios 



53 



Bra and/or Bi"7 fluxes 
H/3 equivalent widths 



9.7 |/m silicate absorption 
UBVRI magnitudes 



JHK magnitudes 



54 



Rosenberg ct al 



Dale ct al. (2007) 
(2008) : Engelbracht et al. (2008) 
Lequeux et al. 



Corbinetal. 

Skillman. Kennicutt. fc Hodge (1989) 7 Lequeux et al. (1979) 
Masegosa, Moles, fc Campos- Aquilar (1994) ; Miller fcllodge (1996) 

van Zee, Havnes. fc Salzer (1997) ; Lee fc Skillman (2004) 
Mendes de Oliveira et al. (2006) : Sidoli. Smith, fc Crowther (2006) 
Rosenberg et al. (2006) : Corbin et al. (2008) 
Engelbracht et al. (2008) ~ 
Begum et al. (2006) : Chengalur et al. (2006) : Paturel et al. (2003) 
Huchtmeier et al. (2005) : Pustilnik fc Martin (2007) 
Hunter fc Gallagher (1985): Hunter, van Woerden. fc Gallagher (1994) 



Lchnert & Hcckman (1995); Annus, Hcckman, & Milcy (1990) 
Gallagher et al. (1991) : Miller fc Hodge (1994) 
Marlowe. Meurer. fc Heckman (1997) : Mendez et al. (1999) 
Gil de Paz. Madore, fc Pevunova*72003) ; James et al. (2004) 
Pustilnik. Pramskii. fc Kniazek (2004) : Hunter fc Elmegreen (20C 



Begum et al. (2006) : Schmitt et al. (2006) 
Rosenberg et al. (2006 , 2008) : Engelbracht et al (2008) 
" French (1980) : Kinman fc Davidson (198T) 
Campbell. Terlevich. fc Melnick (1986 ) 
de Grijp et al. (1992) ; Kennicutt (1992) 
VaccaTTConti (1992) ; Roennback fc Bergvall (1995) 
Miller fc Hodge (1996) ; Izotov et al. (1997) ; van Zee (2000) 
Kong et al (2002): Guseva. Izotov. fc Thuan (2000) 



Kniazev et al. (2003, 2004): Rosenb erg et al. (2006) 
Moustakas fc Kennicutt (2006); Izotov, Thuan, fc Stansihska (2007) 
Moorwood fc Oliva (1988) ; Kawara, Nishida fc Phillips (1989) 
Dale et al. (2001) ; Verma et al. (2003) 
Campbell. Terlevich, fc Melnick (1986) ; Bergvall fc Olofsson (1986) 
Terlevich et al. (1991) ; Kennicutt (1992) ; Vacca fc ContiTl992) 
Roennback fc Bergvall (1995) ; Pustilnik et al. (1999) 
Guseva. Izotov, fc Thuan (2000); van Zee (2000) ; Kong et al. (2002) 



55 

(>3 bands) 



55 

(>1 band) 



Kniazev et al. (2004); Salz er et al. (2005) ; Thuan fc Izotov (2005) 
Izotov et al. (2006) ; Moustakas fc Kennicutt (2006) 
Izotov. Thuan, fc Stansiriska (2007) 
Aitken et al. (1982) ; Phillips. Aitkens, fe Roche (1984) 
Thuan, Sauvage, fc Madden (1999) ; Houck et al. (2004) ; Brandl et al. (2006) 
Salzer. MacAlpine. fc Boroson (1989) : Roennback fc Bergvall (1994) 
Makarova et al. (1998 ); Makarova (1999) 
Mendez fc Esteban (1999) : Cairos et al. (2001) : Bergvall fe Ostlin (2002) 
Paturel et al. (2003) ; Gil de Paz. Madore, fc Pevunova (2003) 
Jangren et al. (2005) : Noeske et al. (2003) 
Pustilnik. Pramskii. fc Kniazek (2004) ; Moustakas fc Kennicutt (2006) 



Spii 



Hunter fc Elmegreen (2006) 
> et al. (1995) ; Vanzi et al. (2000) ; Jarrett et al. (2003) 



Vaduvescu et al. (2005) : Rosenberg et al. (2006) 
Moustakas fc Kennicutt (2006) ; Dale et al. (2007) 
Engelbracht ct al. (2008) 
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